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ABSTRACT 



Aims. We report on the thermal properties and composition of asteroid (2867) Steins derived from an analysis of new Spitzer Space 
Telescope (SST) observations performed in March 2008, in addition to previously published SST observations performed in November 
2005. 

Methods. We consider the three-dimensional shape model and photometric properties derived from OSIRIS images obtained during 
the flyby of the Rosetta spacecraft in September 2008, which we combine with a thermal model to properly interpret the observed 
SST thermal light curve and spectral energy distributions. 

Results. We obtain a thermal inertia in the range 100±50 JK^'m^^s"''^ and a beaming factor (roughness) in the range 0.7-1.0. We 
confirm that the infrared emissivity of Steins is consistent with an enstatite composition. The November 2005 SST theiTnal light 
curve is most reliably interpreted by assuming inhomogeneities in the thermal properties of the surface, with two different regions of 
slightly different roughness, as observed on other small bodies, such as the nucleus of comet 9P/Tempel 1. Our results emphasize that 
the shape model is important to an accurate determination of the thermal inertia and roughness. Finally, we present temperature maps 
of Steins, as seen by Rosetta during its flyby, and discuss the interpretation of the observations performed by the VIRTIS and MIRO 
instruments. 

Key words. Minor planets, asteroids: general - Minor planets, asteroids: individual: 2867 Steins ~ Techniques: spectroscopic 



1. Introduction 

' The Rosetta mission of the European Space Agency flew by as- 
] teroid (2867) Steins on 5 September 2008. Before this flyby, 
Lamy et al. (2008a) and Barucci et al. (2008) performed and 
analyzed Spitzer Space Telescope (hereafter SST) observations 
performed in November 2005, to determine the size, thermal 
properties, and infrared emissivity of the asteroid. At that time, 
the analysis was limited by its unknown size and shape, the un- 
certainties in the photometric properties, and the low signal-to- 
noise ratio (SNR) of the spectra used to derive the emissivity 
[ because its size was much smaller than anticipated. Since the 
, Rosetta flyby, the size and shape (Jorda et al., 201 1), and photo- 
metric properties (Spjuth et al., 201 1) have been determined us- 
ing OSIRIS images. To acquire spectra of higher SNR and better 
. quality, we performed new SST observations of asteroid Steins 
in March 2008, before the Rosetta flyby. With this additional 
data, we can now perform a complete and improved analysis of 
the November 2005 and March 2008 SST data, to determine the 
infrared emissivity of Steins more accurately and refine its ther- 
mal properties (thermal inertia and roughness). These observa- 
tions were not possible with Rosetta since none of its instruments 
cover the spectral range 5-40 yum. 



2. Spitzer Space Telescope observations 

In this article, we use two datasets of SST observations of as- 
teroid Steins as summarized in Table 1. The first dataset was 
obtained on 22 November 2005, when Steins was 2. 13 AU from 



the Sun, and 1.60 AU from SST, at a solar phase angle of 27.2°, 
and is presented in detail in Lamy et al. (2008a). It consists of 14 
spectra of Steins from 5 to 38 fim, taken over one rotation, from 
which Lamy et al. (2008a) derived its thermal light curve and its 
average infrared spectrum. 

The second dataset was obtained on 30 March 2008, when 
Steins was 2.39 AU from the Sun, and 1.88 AU from the SST, 
at a solar phase angle of 23.6°, and is presented for the first time 
in this article. Although Steins was slightly further away from 
the Sun and SST in March 2008 than in November 2005, the 
geometric conditions of the observations are close. We used the 
infrared spectrograph IRS in low-resolution mode to obtain five 
spectra of Steins in the wavelength range 5.2-38.0 fj.m. The short 
wavelength segments (SL2: 5.2-8.5 fim; SLl: 7.4-14.2 yum) were 
acquired with a 60 s ramp, while the long wavelength segments 
(LL2: 14.0-21.5 //m; LLl: 19.5-38.0 //m) were acquired with a 
120 s ramp. We used the red peak-up camera at 22 yum to place 
Steins into the slit with high accuracy. 

The March 2008 data were reduced using standard proce- 
dures developed for IRS, as described in detail in Lamy et al. 
(2008a). The data were first processed by the Spitzer Space 
Center using version S18.0.2 of the data reduction pipeline, to 
produce basic calibrated data (BCD). We then used the BCD to 
extract the spectrum of Steins using version 2.1.2 of the SPICE 
software with standard parameters. According to the SPICE soft- 
ware, the resulting uncertainty in the flux amounts to 3 %. 

After data reduction, a mismatch between the different SL 
and LL modes from 6 to 16 % was clearly visible, depending 
upon the mode. To correct for this elfect, we used the following 
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normalization factors for the different modes: 1 .06 for SL2, 1 .00 
for SLl (our reference), 1.15 for LL2, and 1.20 for LLl. These 
factors are accurate to about 5 %, which increases the overall un- 
certainty in the flux to ~6 %. We believe that the mismatch be- 
tween the different modes, taken at different times, results from 
the changing cross-section of the non-spherical rotating body. 
As we see in section 4, this effect can change the flux by 20 % 
in a few tens of minutes. 

An alternative explanation of the mismatch between the 
SLl and LL2 modes could be the "14 micron teardrop", 
a known effect of the IRS instiTiment, which leads to 
excess emission at the end of the SLl mode (13.2- 
15 yum). However, this effect is difficult to quantify (Spitzer 
Science Center, 2011, IRS Instrument Handbook version 
3.0, http://ssc.spitzer.caltech.edu/irs/irsinstrumenthandbook/). In 
November 2005, the different parts of the spectrum, correspond- 
ing to the different modes (SLl, SL2, LLl, and LL2), were ac- 
quired almost instantaneously relative to the rotation period of 
the asteroid, and could all be smoothly connected without apply- 
ing scaling factors. In particular, no jump was detected between 
SLl and LL2, indicating a negligible "teardrop" in November 
2005. Since we reobserved the same object in March 2008, we 
assume that the "teardrop" is likewise negligible. This assump- 
tion is discussed and justified in Section 5. 

The average spectra of November 2005 and March 2008 are 
displayed in Fig. 1. The flux is lower in March 2008 than in 
November 2005, because of the larger heliocentric and SST- 
centric distances. However, thanks to the longer integration time, 
flie SNR is higher in March 2008. 




15 20 25 
Wavelength {/-im) 



40 



Fig. 1. November 2005 and March 2008 SST spectra of asteroid 
Steins, between 5 and 38 //m. Thanks to the longer integration 
time, the SNR is higher in March 2008. The dashed line rep- 
resents the continuum used to derived the emissivity (Fig. 7), 
which corresponds to a black body at a temperature of 217 K. 



3. Nucleus model 

The interpretation of the infrared data requires a thermal model 
that decribes the energy balance on the surface between the flux 
received from the Sun, the re-radiated flux and the heat conduc- 
tion into the asteroid. This thermal model is associated with a 
shape model that must be used to properly describe the thermal 
infrared light curve. 



3.1. Shape model 

We use the shape model of asteroid Steins constructed by Jorda 
et al. (2011). This model is based on several techniques involv- 
ing non-resolved and resolved images. The portion of the aster- 
oid imaged by Rosetta/ORJSIS is determined by stereo-photo- 
clinometry, while the portion that is not imaged is constrained by 
the inversion of a large set of light curves (Lamy et al. 2008b). 
The overall dimensions of the model are 6.8 x 5.7 x 4.4 km 
as measured along the principal axis of inertia. We use a re- 
sampled version of the model of Jorda et al. (2011), composed 
of 1144 triangular facets. The pole orientation of the model is 
defined by R.A.=90.7 " and DEC=-62.1 °, and the rotation pe- 
riod is 6.04681 ±0.00002 h (Lamy et al. 2008b). Fig. 2 iflus- 
trates the shape model of asteroid Steins as seen from the SST in 
November 2005 and March 2008. 

3.2. Thermal model 

We use the thermal model presented in Lamy et al. (2008a). The 
surface energy balance for each facet of the shape model is given 
by 



(1 - A{X)) — vcos(z) - rjeo-T^ + k— 

rr ox 



(1) 



where A(A) is the Bond albedo as a function of wavelength A, 
FmniX) [Wm^^] is the solar flux as a function of wavelength, 
rh (AU) is the heliocentric distance, v is the illumination factor 
of the facet (v = 1 if the facet is illuminated and v = if the 
facet is in shadow), z is the zenithal angle of the facet, rj is the 
beaming factor introduced by Lebofsky and Spencer (1989), e is 
the infrared emissivity, T (K) is the surface temperature of the 
facet, K (Wm^'K"') is the thermal conductivity, and x measures 
the depth. As the asteroid rotates around its spin axis, the value 
of z changes. 

We solve Eq. (1) using a method similar to that of Spencer et 
al. (1989), described in Groussin et al. (2004). It introduces the 
thermal inertia / - ^KpC, where p (kg/m^) is the bulk density, 
and C (J/kg/K) is the specific heat capacity of the asteroid. We 
use a time step of 7 sec, which is small enough relative to the 
rotation period (~6.05 hr) to ensure relaxation of the numerical 
solution in 100 rotations. As a result, we obtain the temperature 
of each facet as a function of time, over one period of rotation. 
From this temperature profile, we calculate the infrared flux F{?l) 
from each facet as a function of time using Eq. (2) 



-^B(A, T)u cos,{w)dS 



(2) 



where A (km) is the distance to the observer (SST in our case), 
B(A, T) is the Planck function, u is the view factor between the 
facet and the observer (m = 1 if the observer see the facet, m = 
otherwize), w is the angle between the facet and the observer, 
and dS is the surface area of the facet. The total flux received 
by the observer is the sum of the individual fluxes of all facets 
of the shape model. It is computed at different wavelengths to 
construct the spectral energy distribution. 

3.3. Parameters of the thermal model 

Our model has four free parameters: the infrared emissivity e, the 
Bond albedo A(A), the beaming factor //, and the thermal inertia 
/. 

The infrared emissivity e is assumed to be 0.95, the middle 
point of the interval 0.9-1.0 quoted in the literature. Since the 
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Date (start-stop) 


rb 


A 


a 


A\ — A2 


Number 


Integration time 




(AU) 


/"ATA 

(AU) 


(degree) 


Qim) 


of spectra 


(sec) 


22 Nov. 2005 UT 06:23-06:31 


2.13 


1.60 


111 


5.2-38.0 


J 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 07:03-07:10 


2.13 


1.60 


11.1 


5.2-38.0 


J 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 07:38-07:45 


2.13 


1.60 


11.1 


5.2-38.0 


I 


6 (SL mode) or 15 (LL mode) 


22 Nov 2005 UT 08:13-08:20 


2.13 


1.60 


111 


5.2-38.0 


I 


6 (SL mode) or 15 (LL mode) 


22 Nov 2005 UT 08:35-08:42 


2.13 


1.60 


111 


5.2-38.0 


I 


6 (SL mode) or 15 (LL mode) 


22 Nov 2005 UT 09:08-09:15 


2.13 


1.60 


111 


5.2-38.0 


1 


6 (SL mode) or 15 (LL mode) 


22 Nov 2005 UT 09:32-09:39 


2.13 


1.60 


111 


5.2-38.0 


1 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 10:06-10:11 


2.13 


1.60 


111 


5.2-38.0 


1 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 10:38-10:45 


2.13 


1.60 


111 


5.2-38.0 


1 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 11:04-11:12 


2.13 


1.60 


111 


5.2-38.0 


1 


6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 11:35-11:42 


2.13 


1.60 


111 


5.2-38.0 




6 (SL mode) or 15 (LL mode) 


22 Nov 2005 UT 12:12-12:19 


2.13 


1.60 


111 


5.2-38.0 




6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 12:36-12:43 


2.13 


1.61 


111 


5.2-38.0 




6 (SL mode) or 15 (LL mode) 


22 Nov. 2005 UT 13:11-13:18 


2.13 


1.61 


111 


5.2-38.0 




6 (SL mode) or 15 (LL mode) 


30 Mar. 2008 UT 03:05-03:09 


2.39 


1.88 


23.6 


5.2-8.5 (SL2) 


5 


60 


30 Mar. 2008 UT 03:22-03:28 


2.39 


1.88 


23.6 


7.4-14.2 (SLl) 


5 


60 


30 Mar. 2008 UT 03:35-03:47 


2.39 


1.88 


23.6 


14.0-21.5 (LL2) 


5 


120 


30 Mar. 2008 UT 04:00-04:13 


2.39 


1.88 


23.6 


19.5-38.0 (LLl) 


5 


120 



rh: heliocentric distance. A: distance from Spitzer. a: solar phase angle. Ai - A2: spectral range. 



Steins viewed from SST in November 2005 

^ m • m m ^ 

steins viewed from SST m Marcli 2008 



Fig. 2. The shape of asteroid Steins as seen from Spitzer Space Telescope in November 2005 (top) and March 2008 (bottom). Six 
positions are represented, covering one rotation period. In November 2005 the Sun was to the left, and on March 2005 to the right. 



interval is very small and the value near 1 .0, this uncertainty has 
a negligible influence on the calculated thermal flux. 

The Bond albedo A(A) was estimated by Spjuth et al. (201 1) 
using Rosetta OSIRIS images and we use their value. The Bond 
albedo is wavelength dependent and equals to 0.241 at 630 nm. 

The beaming factor 77 follows the strict definition given by 
Lagerros (1998) and therefore only reflects the influence of sur- 
face roughness. In theory, rj ranges from (largest roughness) 
to 1 (flat surface), but in practise must be larger than 0.7 to 
avoid unrealistic roughness, with rm.s. slopes exceeding 45 deg 
(Lagerros, 1998). In this study, rj is unknown and derived from 
the observations. 

For the thermal inertia, we consider six different values 
7=0, 10, 50, 100, 150, and 200 JK-'m-^s'^''^ which cover 
and even extend the range determined by Lamy et al. (2008a) 
of 150±60 JK^'m^^s"'^^. As we demonstrate later (section 4), 
larger values are incompatible with the observations. 



Since e and A(A) are well constrained, we can only vary the 
thermal inertia / and the beaming factor rj. These two parameters 
have a similar efl'ect of changing the surface temperature distri- 
bution. They cannot be determined independently and to each 
value of / corresponds a value for 77. 

4. Analysis of the thermal light curve 

Using the above thermal and shape models, we generate a syn- 
thetic thermal hght curve of asteroid Steins as seen from the SST 
in November 2005. For each value of the thermal inertia 0, 10, 
50, 100, 150, and 200 JK^'m^^s^'^^, we can fit the observations 
to determine the best value for rj. The determination of 77 and the 
phasing of the light curve are performed using a least squares 
technique. The results are presented Table 2. Owing to the 6% 
uncertainty in the flux calibration, 77 is accurate to ±0.05. 

Since we restrict the value of roughness (77) to the range 0.7- 
1.0 as described in section 3.3, the only possible values of the 
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Table 2. Thermal inertia and beaming factor of asteroid Steins 
derived from the November 2005 SST data. 



Thermal inertia / 


Beaming factor 


(JK-'m-^s^'/^) 







1.05 ± 0.05 


10 


1.04 ±0.05 


50 


0.94 ± 0.05 


100 


0.79 + 0.05 


150 


0.69 ± 0.05 


200 


0.61 ± 0.05 



1500 L ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' _ 

; 1=0 JK^'nr^"''^ 7? = !. 05 : 

; 1=10 JK"'m"^s"'/^ 7j = 1.04 : 

; 1 = 50 JK"'m"^s"'/^ 7j = 0.94 I 

; 1=100 JK^'m^^s^'-'^ 7; = 0.79 I 

1400 7 1=150 JK^'m^^s"'/^ r) = 0. 69 Z 




1000 r y, 1 

-0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 

'lime (day) 



Fig. 3. Synthetic thermal light curves of asteroid Steins for dif- 
ferent combinations of thermal inertia / and beaming factor rj 
and the data obtained in November 2005 (diamonds). 

thermal inertia are 0-150 JK~'m~'^s~'^'^, taking into account the 
±0.05 uncertainty in rj. Values of 7=200 JK"'m~^s~^''^ and larger 
are incompatible with the observations as they imply that t] < 
0.7, i.e., an unrealistic roughness. 

Figure 3 illustrates the synthetic light curves of asteroid 
Steins, as seen by the SST on November 2005, for the accept- 
able values of thermal inertia between and 150 JK'^m'^s"'^'^ 
of Table 2. To first order, the agreement is acceptable regardless 
of the thermal inertia, since all the synthetic Ught curves repro- 
duce the observed double peaks, and are consistent with 12 out 
of the 14 data points at the Icr level. However, the first minima 
at ~0.06 days is not well reproduced, the modeled fluxes being 
too low. This first minima corresponds to the part of Steins ob- 
served by OSIRIS and is therefore well defined, so that a shape 
effect seems excluded. We interpret this mismatch between the 
synthetic and observed thermal light curves as resulting from in- 
homogeneities in the thermal properties on the surface. 

Figure 4 illustrates the case where Steins has a thermal iner- 
tia of 100 JK"'m~'^s~^^'^ and two different roughnesses: 77=0.79 
on about 3/4 of its surface and 77=0.72 on the remaining 1/4. The 
transition between these two regions is smoothed out by using 
a combination of the two roughnesses. The second minimum is 




1000 K I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I - 
-0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 
Time (day) 

Fig. 4. Synthetic thermal light curves of asteroid Steins for a 
thermal inertia 7 of 100 JK~^m"^s"^^^ and two different values of 
77, i.e., 0.79 and 0.72. The solid line corresponds to the solution 
that combines two regions with these two different roughnesses. 



well reproduced with 77=0.79, as in Fig. 3, but the first minimum 
is now also well reproduced by instead assuming 77=0.72. The 
variation from 77=0.79 to 77=0.72 between two different regions 
is reasonable as it roughly corresponds to a small change in the 
surface r.m.s. slope of about 10°. These variations have already 
been observed across the nucleus of comet 9P/Tempel 1 for ex- 
ample, there being an average r.m.s. slope of about 12-16° for the 
whole observed surface, but for specific regions slopes of 24-32° 
(Li et al., 2007; Spjuth et al., 2011). 

Similar results are obtained for 7=0 JK^'m^^s"'^^ (respec- 
tively, 10, 50, and 150), where about 3/4 of the surface has 
77=1.05 (respectively, 1.04, 0.94, and 0.69) and the remaining 
1/4 77=0.96 (respectively, 0.95, 0.85, and 0.62). In the latter case 
(7=150), the value of 77=0.62 (±0.05) becomes slightly lower 
than the limit of 0.7 imposed above, but we prefer to be con- 
servative and consider that this solution could also be valid. 

In Fig. 3, each point of the light curve is derived from an 
observed spectral energy distribution (SED), displayed in Fig. 5. 
For the two extreme combinations (1=0, 77=1.05) and (1=150, 
77=0.69) of Table 2, and using our shape model, we compute 
the synthetic SED corresponding to each observed SED. The 
results are illustrated in Fig. 5. Owing to the quasi- similarities 
of the corresponding curves, the two combinations (1, 77) are 
possible and we cannot favor one to the other. We repeat the 
same work for the other intermediate combinations of Table 2 
and find that they are also acceptable. To summarize, the data 
from November 2005 are compatible with a thermal inertia in 
the range 0-150 JK-im-2s-i/2. 

In Lamy et al. (2008a), the shape model was not known with 
sufficient accuracy to perform such a detailed analysis and we 
used the mean of the fourteen individual SEDs (dispayed on 
Fig. 1), with a spherical shape model, thus resulting in slightly 
different results (1=150±60 JK^'m^^s"'^^). This emphasizes the 
importance of using an accurate shape model to interpret in- 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 
Time (day) 

Fig. 6. Synthetic thermal hght curves of asteroid Steins for 
different combinations of thermal inertia / and beaming fac- 
tor T] and the data obtained in March 2008 (diamonds). Each 
data point corresponds to the March 2008 integrated spectrum 
displayed in Fig. 1, normalized by the factors introduced in 
Section 2: 1.06 for SL2, 1.00 for SLl (our reference), 1.15 for 
LL2, and 1.20 for LLl. 



frared thermal light curves to properly differentiate the effects 
of the shape from that of the thermal parameters. 

Once we have determined the thermal parameters for 
November 2005, we can check their consistency with the March 
2008 data. Figure 6 illustrates the synthetic thermal light curve 
of asteroid Steins as seen from the SST for the valid combina- 
tions (0 < / < 150 JK-im-2s-i^2) derived above (Table 2). The 
spectra from March 2008 were taken at different times in differ- 
ent modes (see Table 1), thus providing information on the flux 
variations due to the rotation of the asteroid. Its rotation period 
is not accurate enough to obtain the absolute phase relative to 
November 2005 and the phase shift was adjusted to match the 
March 2008 data. According to the normalization factors pre- 
sented in section 2 for the different modes and as illustrated in 
Fig. 6, we most likely observed the maximum cross-section. The 
low thermal inertia of 0-10 JK~'m~^s~'^^ are unambiguously ex- 
cluded. Therefore, the thermal inertia of Steins must lie in the 
range 50-150 TK-^m-h-^'^, that is 100±50 TK-^m-h-^'^, in 
agreement with the March 2008 data. 

5. Infrared emissivity 

The mid-infrared spectral domains contain features that are diag- 
nostic of the surface composition, such as the Christiansen peaks 
(8-9.5 i-im), Reststrahlen bands (9-12 fim and 14-25 fim), and 
Transparency features (11-13 /im). Using the November 2005 
SST data, Barucci et al. (2008) detected these features on as- 
teroid Steins, and by comparing the observed spectra with lab- 
oratory measurements, concluded that the thermal emissivity 
of Steins is similar to those of enstatite achondrite meteorites, 
which is consistent with Steins being an E-type asteroid. 



C R T 




CRT 

0.2 ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ - 

5 10 15 20 25 30 35 40 

Wavelength (/^m) 

Fig. 7. Spectral variation in the infrared emissivity of asteroid 
Steins derived from the March 2008 observations, compared to 
that derived from the November 2005 observations and to two 
laboratory measurements for enstatite and aubrite (Barucci et 
al., 2008). In the case of the March 2008 data, two spectra are 
represented, the nominal spectrum and an alternative spectrum 
incorporating a "teardrop" effect (see text). 



Figure 7 illustrates the infrared emissivity of asteroid Steins 
derived from the March 2008 data, which is obtained by divid- 
ing the observed spectra by the continuum. In our case, the con- 
tinuum is a black body with a color temperature of 217 K, as 
displayed in Fig. 1. The infrared emissivity of asteroid Steins 
obtained in March 2008 is compared to that of November 2005 
and to laboratory measurements for enstatite and aubrite. The 
Christiansen (C), Reststrahlen (R), and Transparency (T) fea- 
tures are clearly visible in March 2008. The correspondence be- 
tween both the C and R peaks and those of enstatite and aubrite 
is much better than for the November 2005 case, because of 
the higher SNR. Moreover, the broad emission at ~15 /im is 
also more accurately reproduced. It reinforces the conclusion of 
Barucci et al. (2008) that Steins has an enstatite composition. 

As explained in Section 2, we assume a negligible "14 mi- 
cron teardrop" for the nominal spectrum of March 2008. To 
justify this assumption, we display in Fig. 7 the same spec- 
trum, but assuming a non-negligible "teardrop", i.e., ignoring 
the SLl flux beyond 13.2 jum. In this case, the Christiansen and 
Reststrahlen peaks are still present, but the Transparency fea- 
ture disappears, which contradicts the results of November 2005, 
when the "teardrop" was negligible. The consistency between 
the November 2005 and March 2008 emission spectra justifies 
our assumption. 
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6. Temperature maps of asteroid Steins at the time 
of the Rosetta f lyby 

We now consider the temperature maps of Steins as predicted 
by our model for the 5 September 2009 Rosetta flyby to discuss 
the interpretation of the observations performed by the VIRTIS 
(Coradini et al., 1999; Coradini et al., 2008) and MIRO (Guilds 
et al., 2007; Guilds et al., 2010) instruments. We perform the 
computation at two different times, 10 minutes before and at 
closest approach, and for two different combinations of thermal 
inertia and roughness: (i) 7=0 JK 'm^'s^'^^ and ?7=1.05, and (ii) 
7=100 JK-^m-^s"^/^ and t]=0J9. Figure 8 illustrates the results. 
These maps are available in digital form on request from the au- 
thors. Steins rotates by about 10° in 10 minutes inducing only a 
sUght change in the surface temperature distribution. However, 
the viewing angles are quite different, revealing different regions 
of the surface. 

Any differences between the two cases (i) and (ii) is not re- 
flected in the maximum temperature, 252 K for JK^'m^^s^'^^ 
and 77=1.05 and 245 K for 100 JK-^rirh-^/^ and 77=0.79, but 
prominently in the temperature distributions. In the case of a null 
thermal inertia, there is instantaneous adaptation of the surface to 
the solar insulation and the temperature drops to K in the shad- 
owed regions and night side (Fig. 8, middle panel). This leads to 
large temperature gradients on the surface close to the termina- 
tors, and a maximum day/night amplitude. For a thermal inertia 
of 100 JK^'m^^s"''^, the amplitude of the day/night temperature 
variations is smaller than for 0, and the night side temperature 
only drops to about 150 K, close to the equator. Moreover, the 
adapation to solar insulation is not instantaneous as for a null 
thermal inertia, the highest temperature location being shifted to 
the afternoon and the temperature rising less rapidly in the morn- 
ing hemisphere (Fig. 8, bottom panel). These are well known ef- 
fects that, combined with temperature measurements of the sur- 
face obtained by VIRTIS or MIRO, can more tightly constrain 
the thermal properties of Steins. 

The MIRO instrument observed Steins in the millimeter 
(1.6 mm) and submiUimeter (0.53 mm) wavelength ranges 
(Guilds et al., 2010). From their measurements, which are dif- 
ficult to interpret because of pointing uncertainties, Gulkis et al. 
(2010) derived a constraint on the average emissivity of 0.85-0.9 
at 1 .6 mm and 0.6-0.7 at 0.53 mm. They also give a range of 450- 
850 IK'^m'^s"^**^ for the thermal inertia but, as conceded by the 
authors, MIRO itself provides very little direct constraints on the 
thermal inertia, and this range is only based on their interpreta- 
tion of the VIRTIS (unpublished) observations. As confirmed by 
a discussion with the authors, the MIRO measurements alone 
carmot rule out a lower thermal inertia of 100-150 JK"^m"^s"'^^ 
as we found. 

The VIRTIS observations give a maximum temperature of 
about 230 K close to the sub-solar point (Coradini et al., 2008), 
implying a thermal inertia of 140 JK"^m"^s"''^ with 77=1.0, 
consistent with our results. According to Fig. 6 of Gulkis et 
al. (2010) showing unpublished VIRTIS temperature measure- 
ments at different local solar phase angles, the diurnal temper- 
ature variation is too small to be compatible with such a ther- 
mal inertia, and a larger value is required. In this context, we 
note that a similar problem was encountered for the nucleus 
of comet 9P/Tempel 1 (Groussin et al., 2007) and it was con- 
cluded that this effect most Ukely results from a small-scale sur- 
face roughness (not taken into account when generating Fig. 8), 
which tends to increase the flux at larger local solar phase an- 
gles compared to a smooth surface (Rozitis and Green, 2010). 
This interpretation clearly differs from that proposed by Gulkis 



et al. (2010), but until a full and comprehensive analysis of the 
VIRTIS observations of Steins is performed and published, we 
cannot solve this issue. 

7. Conclusions 

We acquired new thermal infrared spectra of asteroid Steins with 
the SST in March 2008 and have interpreted these data, in addi- 
tion to data obtained in November 2005, using the shape model 
and photometric properties inferred from the data of the Rosetta 
flyby of 5 September 2008. The two datasets are consistent with 
a thermal inertia of 100±50 JK^'m^^s"'^^ and a beaming factor 
(roughness) of 0.7-1.0 (Table 2). This range of thermal inertia 
is barely compatible with the constraint 7 <100 JK~^m~^s"''^ 
obtained by Delbo et al. (2007) from their compilation of main- 
belt asteroids (MBAs) larger than 100 km, but consistent with 
the results 7 <400 JK^'m^^s^'^^ (Nesvomy & Bottke, 2004) and 
7 <200 JK-im-2s-i/2 (Harris et al., 2009) for the Karin fam- 
ily, whose mean size of about 4 km is comparable to that of 
Steins. On the basis of the thermal properties of 17 MBAs in the 
size range 25 -1000 km, Delbo et al. (2007) and Delbo & Tanga 
(2009) suggested that the thermal inertia increases with decreas- 
ing size. Our results for Steins do not confirm this putative corre- 
lation for kilometer- size MBAs as it would have implied a ther- 
mal inertia of several thousands JK 'm^^s for Steins based 
on the correlation law of Delbo & Tanga (2009). 

We confirm that the infrared emissivity of Steins is consistent 
with an enstatite composition, in agreement with the results of 
Barucci et al. (2008). 

The November 2005 SST thermal light curve is more accu- 
rately interpreted by assuming inhomogeneities in the thermal 
properties of the surface, or more precisely two different regions 
with different roughness. A smaU change in the beaming factor 
of about 0. 1 is required between the two different regions, which 
is compatible with observations of other small bodies, for exam- 
ple the nucleus of comet 9P/Tempel 1 (Li et al., 2007; Spjuth et 
al., 2011). 

Overall, our results emphasize that the shape model is im- 
portant to an accurate determination of the thermal inertia and 
roughness. 

However, uncertainties remain, in particular for the spatial 
variations in the thermal properties, which requires in-situ mea- 
surements of the surface temperature distribution, as illustrated 
for the nucleus of comet 9P/Tempel 1 (Groussin et al. 2007). We 
hope to make progress in this direction with the flybys of aster- 
oid 21 Lutetia in July 2010 by Rosetta, of comet 103P/Hartley 2 
in November 2010 by EPOXI, and even more promising with the 
Rosetta rendez-vous with comet 67P/Churyumov-Gerasimenko 
in 2014. 
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Fig. 5. Individual spectral energy distribution (SED) observed in November 2005, with their corresponding synthetic SED for the 
extreme combinations of thermal parameters : blue line 1=0 JK^'m^^s"'^^, 77=1.05; red line 1=150 JK^'m^^s"'^^, 77=0.69. 
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Fig. 8. Asteroid Steins as seen from the Rosetta spacecraft during the 5 September 2008 flyby. The left column corresponds to 10 
minutes before closest approach, while the right column corresponds to closest approach. The top panel displays images acquired 
by the OSIRIS camera in the visible. The other two panels display temperature maps: the middle panel corresponds to the case 
7=0 JK-im-2s-i/2 and //^l.OS, and the bottom panel to 7=100 JK-'m-^g-i/^ and ?7=0.79. 



